Dynamic projection on Feshbach molecules: a probe of pairing and phase fluctuations 
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We describe and justify a simple model for the dynamics associated with rapid sweeps across a 
Feshbach resonance, from the atomic to the molecular side, in an ultra cold Fermi system. The 
model allows us to relate the observed molecule momentum distribution, including its dependence 
on the sweep rate, to equilibrium properties of the initial state. For initial state near resonance, 
we find that phase fluctuations sharply reduce the observed condensate fraction. Moreover, for 
very fast sweeps and low temperatures, we predict a surprising nonmonotonic dependence of the 
molecule condensate fraction on detuning, that is a direct signature of quantum phase fluctuations. 
The dependence of the total molecule number on sweep rate is found to be a sensitive probe of 
pairing in the initial state, whether condensed or not. Hence it can be utilized to establish the 
presence of a phase fluctuation induced 'pseudogap' phase in these systems. 



Experiments with ultra cold fermions near a Feshbach 
resonance (FR), opened a new window to the study of 
strongly interacting condensates^, 0, 0- They can ac- 
cess the strongly coupled regime intermediate between 
weak pairing BCS superfluidity and a BEC of molecules. 
In this region, deviations from mean field theory a la BCS 
El 00] are expected to be large. In addition, the abil- 
ity to rapidly vary the interaction parameters, provides a 
unique opportunity to study quantum dynamics far from 
equilibrium|3,|a,l2|- The focus of recent experiments has 
been to utilize dynamics, in this case rapid magnetic field 
sweeps across the FR, in order to probe the equilibrium 
properties of the condensate in the crossover regionj^Q- 
The idea was to convert cooper pairs in the initial state, 
which would otherwise unbind during free expansion, into 
molecules. The fact that the molecule momentum distri- 
bution depended on the start position has been offered 
as strong indication that indeed properties of the initial 
equilibrium state were being accessed. However, precise 
connections to such properties are lacking. 

In this letter we formulate an approach that allows us 
to relate equilibrium properties of such strongly inter- 
acting fermions to the measurements in the dynamical 
experiments. Our approach targets the regime of 'fast' 
magnetic field sweeps. Roughly, these are ramp rates at 
which the conversion efficiency of atoms into molecules 
changes significantly with changing ramp rate. Empiri- 
cally in the K system this implies a ramp rate faster 
than 1/20 G//zsec|3, while it was too fast to be accessed 
in 6 Li|3|]. The physics in this regime turns out to have a 
remarkable simplicity, and for this reason we will focus 
on it in this letter, although the bulk of experiments to 
date have been done at slower ramp speeds. It is hoped 
that the results of the present study will stimulate exper- 
imental investigation of this regime. 

Consider for a moment the extreme limit of an in- 
finitely fast sweep. Then, we can apply the sudden ap- 
proximation and simply use the initial state to evaluate 
the final population of molecules. This intuitive picture 
was advocated by Regal et a/Q (see also 0). Using 
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FIG. 1: Model dynamics. Field sweep is effectively sudden up 
to detuning i/* depending on the sweep rate v. At this point 
the scattering length (and molecule size) is a*. The molecule 
is assumed to evolve adiabatically from there. 

the same assumption, Diener and Ho 01 estimated the 
fraction of condensed molecules as a function of the de- 
tuning of the initial state from resonance. The molecule 
momentum distribution in this approach is given by 
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where (fif(k) is the molecular WF at the end of the ramp, 
deep in the molecular side of the resonance. The aver- 
age ( ) i is taken over the initial state, which in Ref. 
was assumed to be a BCS WF, calculated within the 
mean field theory of 0. Indeed, this gives a bimodal 
momentum distribution, with a peak at q = due to 
cooper pairs projected to molecules; as well as molecules 
at q 7^ due to pairwise projection of atoms belonging to 
different cooper pairs. Although this is roughly what is 
seen in the experiments, it is unsatisfactory in two impor- 
tant respects, (i) Mean field theory is used to calculate 
the initial state, whereas it is well known that in the 
crossover regime near the resonance, quantum and ther- 
mally induced phase fluctuations play a very important 
role, (ii) The approach ignores the dynamical aspects of 
the experiment and is by definition unable to predict the 
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dependence of measured quantities (eg. conversion effi- 
ciency and molecular condensate fraction) on ramp rate. 

Both these issues are addressed in this paper. The first, 
by going beyond the mean field theory, using the RPA 
[l3| to include Gaussian phase fluctuations (or "non con- 
densed cooper pairs"). In the BCS to BEC crossover 
regime these sharply reduce the observed condensate 
fraction relative to the mean field result of Rcf. 
even at zero temperature. For fast sweeps at low tem- 
perature this leads us to predict a surprising non mono- 
tonic behavior of the molecule condensate fraction versus 
detuning of the initial state from resonance. Point (ii) 
raised above is addressed by an effective model for the 
dynamics. For fast sweeps we argue that the time evo- 
lution of the system can be approximated by a 'sudden' 
part followed by an 'adiabatic' time evolution part. The 
point at which the time evolution changes character de- 
pends on the ramp rate, the sudden evolution persisting 
for longer at higher ramp rates. Thus, the dynamics is 
approximated by projecting the initial state of the sys- 
tem onto an effective molecular WF, which is determined 
by the ramp rate. So, we will be using Eq. l[T]l. but with 
the effective molecular WF, ip that depends on the ramp 
rate. In this way we are able to obtain the paramet- 
ric dependence of various measured quantities, such as 
conversion efficiency into molecules and molecular con- 
densate fraction on the ramp rate, and the parameters 
of the initial state. At these fast sweep rates, we find 
that conversion of Cooper pairs into molecules is vastly 
more efficient than that of uncorrelated pairs of atoms. 
This applies both to condensed and uncondensed Cooper 
pairs, and ultimately results from a short distance sin- 
gularity present in the Cooper pair wavefunction (WF), 
that allows them to have a non-negligible overlap even 
with small sized molecules. Since these dynamical mea- 
surements are equally sensitive to noncondensed Cooper 
pairs, they can potentially probe the phase fluctuation 
induced pseudogap phase (where thermal and quantum 
phase fluctuations destroy the condensate but pairing re- 
mains). Details such as the momentum distribution of 
noncondensed pairs may also be accessed. 

We now motivate the model dynamics via physical ar- 
guments and derive the consequences for the system of 
interest. Later we perform a nontrivial check by showing 
that this simple scheme indeed reproduces the physics in 
a nontrivial toy model (the Dicke model) which is solved 
numerically without approximations. In what follows we 
concentrate on a wide FR, relevant to experiments in 
4o i4T0 and 7 Li^ ]. That is, g s = g^j2 » e f where g 
is the coupling between the open and closed channels and 
n is the atom density. For most purposes it is then possi- 
ble to neglect the occupation of closed channel molecules. 
The problem reduces to spin-^ fermions interacting via 
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Ultraviolet divergences are avoided in the standard fash- 
ion by writing all physical results in terms of an s-wave 
scattering length a instead of u. Near resonance the scat- 
tering length diverges as a w — mg 2 /Amy where m is the 
atomic mass, and v is the detuning in energy units, which 
is related to the magnetic detuning via the magnetic mo- 
ment difference A/j, between closed and open channels: 
v = Afi(B - B ). 

Our two stage approximation of the dynamics is de- 
picted graphically in Fig. ^ For slow sweeps, atoms 
are converted adiabatically to weakly bound Feshbach 
molecules (FM) The binding energy of these 

molecules is Eb = —l/(ma 2 ) oc v 2 (for kfd << I). Now 
consider a change in the detuning parameter at a con- 
stant rate v. Once the binding energy is large enough, 
such that Ei, <C E 2 , the time evolution is expected to be 
adiabatic. So, for a given sweep rate u, there is a char- 
acteristic detuning i/*, which marks a dynamic crossover 
for the system. At i/*, E^ w E 2 ^^l- O ur simplified two 
stage model for the dynamics approximates the change 
up to as sudden, while the rest is considered as per- 
fectly adiabatic. Thus the WF at the initial state is ef- 
fectively "projected" on FMs that occur at detuning v±. 
These molecules evolve adiabatically into more tightly 
bound ones corresponding to the final value of the field, 
while the rest of the atomic population remains unbound. 

Within this model, a faster sweep rate simply moves 
the break point v+ to larger negative detuning, leading to 
effective projection of the initial WF on smaller FMs. We 
can establish a precise connection between the sweep rate 
and the size of the effective molecule WF. The typical size 
of the FM at v+ is the scattering length a*(i^*). Using 
Ef, = — I /(ma 2 ), and the relation between the scattering 
length and the detuning v quoted earlier, with the break 
condition E 2 (a±) — dEj,(a i ,) / dt we obtain 
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We now employ the two stage model and Eq. to 
determine the main features of the final molecule distri- 
bution and its dependence on both initial state and sweep 
rate. We need to evaluate the following correlation func- 
tion, in the initial state: 



n m (q) = I rf 3 ^ 3 ^V*(fc)<p(fc')(4+kT4-ki c f-k'i c f+k'T)- 

(4) 

The size of the molecular pair WF (p(k) is a+, as pre- 
scribed by the model dynamics. First consider a mean 
field (BCS) approximation of the initial state. In this 
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case the correlation function appearing in (J2J may be 
factorized to obtain, as in Ref. [ill: 
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The first contribution, proportional to the anomalous ex- 
pectation value, gives the condensed part of the distribu- 
tion. In fact this is just the square of the overlap between 
the final molecule WF and the "cooper pair WF". The 
second, non condensed part, contains the normal expec- 
tation values. To make further progress we note that 
if the molecule size a* is much smaller than the inter- 
particle spacing, we can replace the exact WF with a 
box WF in momentum space, of the same spatial extent. 
We take tp{k) = v /3/4 7 ra 3/2 for k < 1/a* and <p(k) = 
outside this sphere. Now (p(k) serves as a cutoff to the 
relative momentum integrals. Using this in JJjJ, the num- 
ber of normal molecules is found to be 



N n = {N/2)(k f a+f 
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where N is the total atom number. This result is easily 
understood in terms of the overlap of a random pair of 
atoms in the Fermi sea, with a molecular WF. The con- 
version efficiency is then proportional to the ratio of the 
molecule volume to that occupied on average by nearby 
atoms. In contrast, the number of condensed molecules 
calculated from J5j is 
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The last step relies on the high momentum divergence of 
the integral, due to which it depends crucially on the high 
momentum cutoff. The result (0 will have important 
consequences on interpretation of experiments, and it is 
worthwhile to remark on its physical origin. The number 
of condensed molecules is proportional to the square of 
the overlap of the molecular WF ip m with the Cooper pair 
WF ip c . Now, the molecule WF is appreciable only within 
a region r < a*, where ip m (r) ~ a* 3 ^ 2 . Since this region 
is smaller than the average inter-particle distance, the 
molecules probe the Cooper pairs at very short distances. 
Now, the Cooper pair WF has a singular short distance 
behavior, ¥? c (r) ocl/r. Therefore \ (ip m \ip c )\ 2 oc a* which 
is a much larger overlap than might be naively expected. 
Note, in solid state systems where a natural short dis- 
tance cutoff exists, these features are absent and hence 
have not been emphasized. 

Combining the mean field results for Nq and N n we 
can evaluate the condensate fraction in the molecule dis- 
tribution: 
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The dependence on the initial state enters this expres- 
sion through the factor e//A. For large positive detun- 
ing we expect weak pairing A/ej ~ e~ l ' k f a , while close 
to resonance A ~ e*. For sufficiently fast sweep rates, 
kfa* << 1 we have /mf ~ 1 in the vicinity of the res- 
onance. This is because Cooper pairs are much more 
efficiently converted into molecules, and within the mean 
field approximation used above, all Cooper pairs are con- 
densed. However, especially in this region close to the 
resonance we expect phase fluctuations (Cooper pairs at 
finite momenta) to be excited, which will lead to / < 1. 

To obtain such corrections to the mean field result (0 
we calculate the correlation function Q within the RPA 
approximation. The details of the calculation are left to 
|l2j, here we briefly outline the main steps and the re- 
sults. Following Ref. we consider a path integral 
representation of the partition function Z determined by 
the Hamiltonian [5] The interaction term may be decou- 
pled with a Hubbard-Stratanovitch pair field A(q, lu). In 
order to extract the desired correlation functions, we in- 



troduce a source term J(q, uS) f d 4 kip(k)c^ Q ^ 

h.c, where we have used the four vector notation Q = 
(q, u)). At T > the integral over imaginary frequencies 
is converted into the usual Matsubara sum. Then, the 
desired molecular distribution function is: 
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While the saddle point approximation gives us the 
BCS result, here we go beyond and expand A(q, u>) — 
Ao^(Q) + 7?(q, w), where A is the saddle-point value 
of the gap. The RPA approximation involves integrating 
out the fermions and expanding the resulting effective ac- 
tion to quadratic order in the rj and J fields. The poles 
in the r\ propagator give the collective mode spectrum. 
Finally, integrating out the r\ fields gives us an effective 
action solely in terms of J from which the molecular dis- 
tribution function n m (q)MF + 5n rn (q) (Sn(q) is the con- 
tribution due to fluctuations) may be evaluated. 

Assuming that the fluctuation contribution is domi- 
nantly from superfluid phonons, we can make a small mo- 
mentum, small frequency expansion, where the J propa- 
gator takes the form F(fc/a Jr , kfa)[(cq) 2 — (iiv) 2 }^ 1 , where 
c is the sound velocity in the superfluid and the func- 
tion F will be discussed shortly. It is important to note 
that the implied linear dispersion is only an approxima- 
tion due to the low q and lo expansion. At momenta 
q ^ 4mc the dispersion becomes quadratic. In addition 
at u> > A phase fluctuations can decay into quasiparti- 
cle excitations, leading to damping in the BCS limit at 
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Carrying out the Matsubara summation we obtain the 
RPA correction to the molecule momentum distribution 

r / n /, , ,coth(co/2T) , s 

5n(q) = F(fc/a*, k f a) '-. (10) 
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FIG. 2: Calculated molecule condensate fraction versus 
— 1/fe/a in the initial state (bottom axis),and versus mag- 
netic field detuning for parameters of 40 K 2] (top axis). The 
sweep rate corresponds to fe/a* = 0.3. Solid lines include the 
effect of quantum and thermal phase fluctuations. Dashed 
line shows the (T=0) mean field result for comparison. 



The number of non-condensed molecules due to these 
fluctuations is found by integrating over q up to a natu- 
ral cutoff. As discussed above, such a cutoff is provided 
in the BCS limit by q% = l/£, and in the BEC limit by 
(72 = 4mc. To cover the whole range we use the cut- 
off q^ 1 = Qi 1 + q^ 1 which interpolates between the two 
limits. 

The variation of the above fluctuation correction with 
sweep rate is encoded into the dependence of the func- 
tion F{kfa i , 1 kfa) on fc/a*. The leading dependence on 
fc/a* is found to be linear, which arises from the fact 
that non-condensed Cooper pairs are as efficiently trans- 
formed into molecules as condensed Cooper pairs. By 
contrast, conversion of unpaired atoms into molecules is 
much less efficient and scales as (fc/a*) 3 . 

The full function F(kfa+, kfO) can be computed along 
the crossover from BCS to BEC|f2|. Fig. El depicts the 
calculated condensed molecule fraction including fluctu- 
ation corrections for sweep rate corresponding to projec- 
tion at kfd-i, — 0.3. The nonmonotonic dependence on 
detuning toward the BCS limit is due to competition be- 
tween two effects. On the one hand decreasing phase 
fluctuations act to increase the condensed molecule frac- 
tion. On the other hand the ratio e//A appearing in 
the normal atom contribution diverges exponentially at 
large detuning and eventually leads to vanishing of the 
molecular condensate fraction. 

We now turn to a non trivial check of the two stage 
approximation of the dynamics. Based on this model we 
found enhanced conversion efficiency at rapid sweeps of 
cooper pairs (both condensed and non condensed) into 
FMs, compared to unpaired atoms. Using 0, and the 
linear scaling in /c/a* of the number of molecules arising 
from cooper pairs, we obtain a parametric dependence 
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FIG. 3: Simulated mean Dicke model dynamics: Final 
molecule number as a function of the sweep rate, (a) Initial 
state close to resonance (strong pairing). The curves for dif- 
ferent coupling constants g s fall on a universal curve oc Xjv 1 ^ 3 
at fast rates, in agreement with and 0. (b) Initial state 
far detuned (v = 400e/, g s = 6e/), i.e. weak pairing. The 
1/v behavior due to @ dominates. 

on the sweep rate oc (gl/v) 1 ^. By contrast, from © 
and J2J), the conversion efficiency for unpaired atoms is 
much lower, oc gl/v. We now wish to verify these non- 
trivial dependences using a numerical simulation of the 
dynamics, not relying on the two-stage model. We start 
from the microscopic two channel HamiltonianQ 

H = y^(ek - M)c* k c gk + y^(e fc /2 + v - 2^)m\m k 

kcr k 
kq 

which is regularized by absorbing the high momentum 
cutoff into renormalized detuning parameter [(j . Note, 
TO q describes a (bare) closed channel molecule. 

We compute the dynamics of HII|) within the Dickie 
model (i.e. keeping only m q= o). The initial state is taken 
to be the equilibrium solution at detuning v > 0. Then 
v is changed at a constant rate to far negative detuning 
where the equilibrium population of mo would be ~ 96% 
and the resulting dynamics (see e.g. 0) is calculated 
numerically with no approximation. Following the ex- 
periments we count the number of molecules in the final 
state, which to a good approximation is |(mo)| 2 . Accord- 
ing to the two stage model this is directly related to the 
number of FMs produced after the sudden stage. 

The results are summarized in Fig. |3| Most impor- 
tantly, for initial state at resonance (i.e. strong pairing), 
the dependence on sweep rate fits 1/v 1 / 3 , as expected 
from our two stage model. By contrast, for weak pair- 
ing, deep in the BCS regime the normal contribution l/z> 
dominates (Fig. ^p). The strong suppression of the pref- 
actor of 1/v is an artifact of the Dickie model, which 
prohibits occupation of molecules with q =/= 0. 

In summary, we presented a simple model of the dy- 
namics of rapid magnetic field sweeps across the FR. This 
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allowed to relate the measurements of such dynamical ex- 
periments to properties of the initial equilibrium state of 
the Fermions. For rapid sweeps at low temperatures in 
the crossover region we predict a a non monotonic behav- 
ior of the final condensed molecule fraction versus detun- 
ing, which is a direct signature of the quantum phase 
fluctuations in the initial state. In addition, the conver- 
sion efficiency to molecules at fast sweep rates is argued 
to be a sensitive probe of pairing, whether in a condensed 
state or not. This can be used to establish the presence 
of a fluctuation induced pseudogap phase. In this regime 
we expect enhanced conversion efficiency due to the pres- 
ence of noncondensed fermion pairs, despite a vanishing 
condensate fraction. 
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